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Abstract 

The model under consideration is a two-dimensional two-component plasma, 
i.e., a continuous system of two species of pointlike particles of opposite charges ±1, 
interacting through the logarithmic Coulomb interaction. Using the exact results for 
the form-factors of an equivalent Euclidean sine-Gordon theory, we derive the large- 
distance behavior of the pair correlation functions between charged particles. This 
asymptotic behavior is checked on a few lower orders of its /3-expansion (/3 is the 
inverse temperature) around the Debye-Hiickel limit /? — > 0, and at the free-fermion 
point /3 = 2 at which the collapse of positive-negative pairs of charges occurs. 
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1 Introduction 



We consider a two-dimensional (2D) two-component plasma (TCP), i.e., a neutral con- 
tinuous system of two species of pointlike particles of opposite charges ±1, interacting 
through the 2D logarithmic Coulomb interaction. Classical equilibrium statistical me- 
chanics is used. The system is stable against collapse of positive-negative pairs of charges 
for the dimensionless coupling constant (inverse temperature) (3 < 2. In the region /3 > 2, 
one has to attach to particles a hard core in order to prevent the collapse. In the limit 
of a small, but nonzero, hard core, the Kosterlitz-Thouless phase transition from a high- 
temperature conducting phase to a low-temperature insulating phase [0] takes place at 
around /3 = 4. 

The 2D TCP is the first continuous fluid in dimension higher than one with an exactly 
solvable thermodynamics, in the bulk as well as at the surface for specific boundary 
conditions. The bulk equation of state has been known for a very long time 0]. The 
other bulk thermodynamic properties (free energy, internal energy, specific heat, etc.) 
have been obtained in the whole stability range /3 < 2 of the plasma in a recent paper 
0. The mapping onto a bulk 2D Euclidean (classical) sine-Gordon field theory with a 
specific (conformal) normalization of the cos-field was made, and recent results about 
that field theory [|] - @, were used. In subsequent works, the surface tension of the same 
model in contact with an ideal conductor and an ideal dielectric rectilinear walls 
was obtained. A mapping onto a boundary 2D sine-Gordon field theory with a Dirichlet 
and a Neumann boundary condition, respectively, was made, and known results |^ about 
those integrable boundary field theories [|lOl were applied. 



In this paper, we derive the large-distance behavior of the bulk particle correlation 
functions for the 2D TCP by exploring the form-factor theory of the equivalent sine- 
Gordon model. 

Exact sum rules are known for the correlation functions. The zeroth and second mo- 
ments of the charge correlation function have the well-known Stillinger-Lovett values |TT 



The zeroth moment of the density correlation function is related to the compressibil- 



ity (here exactly known), and the second moment has been recently obtained [T^, [13 
An exact analysis of pair correlations in the 2D TCP can be done in two cases: in the 
high-temperature Debye-Hiickel limit /3 — > and just at the collapse point (3 = 2 which 
corresponds to the free-fermion point of an equivalent Thirring model (although, at a given 
fugacity, the free energy diverges, truncated particle distributions are finite at /? = 2 and, 
as is generally believed, also in the whole range /3 < 4). We have not found in the lit- 
erature an exhaustive treatment of the limit /5 — > 0, so we study this limit in detail in 
the present work by using the technique of a renormalized Mayer expansion in density 
|1I3| , p^ . At /3 = 2, the exact forms of the bulk truncated distributions (Ursell functions) 
were found in ref. via a continualization of the Gaudin's model |16|. The Ursell 
functions were computed at /3 = 2 also for inhomogeneous situations when the 2D TCP, 
being confined to domains of various geometry, is in contact with a charged hard wall. 



polarizable interface [T^], an ideal conductor wall ||T^ - |1T9[ and an ideal dielectric wall 
|]20| , A short-distance expansion of pair correlations can be done in principle for an 



arbitrary (3 by using Dotsenko's integrals |22]; interestingly, correlations between equally 



charged particles change their analytic structure at short distance at /5 = 1 [23 



In the 2D sine-Gordon theory, a differential equation for generalized correlation func- 
tions was derived at the free-fermion point [Q. The classical 2D sine-Gordon model 
can be regarded as a conformal (namely Gaussian) field theory, perturbed by the cos- 
field. For such theories, the short- distance expansion for multipoint correlation functions 
can be systematically obtained by using the Operator product expansion [^, combined 
with the Conformal perturbation theory [^]. The most efficient way to study the large- 
distance behavior of correlations is provided by the form-factor approach. For the 2D 
sine-Gordon model, the form factors of some local operators for all kinds of particles the 
theory contains were evaluated in a series of papers - |3^. In particular, the form 
factors for the lightest particle (elementary breather), dominant when this particle exists, 
are analytic continuations of their counterparts in the sinh-Gordon theory [^, 0. In the 
sinh-Gordon model, the corresponding particle is the only massive one; the contributions 
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of all form factors were summed up and a closed expression for any correlation function 



as a determinant of an integral operator was obtained in ref. [33 



The paper is organized as follows. In section 2, we introduce the notation and briefly 
review some important information on the mapping of the 2D TCP onto the 2D sine- 
Gordon model 0. Section 3 deals with the form-factor theory for the sine-Gordon model, 
written in a way accessible to non-specialists as well. The large- distance asymptotic of 
pair particle distributions for the 2D TCP is presented. The obtained results are checked 
around the Debye-Hiickel limit f3 and at the free fermion (3 = 2 point in section 4. 
A brief recapitulation and some concluding remarks are given in section 5. 



2 Sine-Gordon representation of the plasma 

We consider an infinite 2D space of points r G i?^. The 2D TCP, defined in this space, 
consists of point particles {j} of charge {cTj = ±1}, immersed in a homogeneous medium of 
dielectric constant = 1. The interaction energy of particles is given by Y.i<j o'iO'jv(\ri — rj\), 
where the Coulomb potential v is the solution of the 2D Poisson equation 

Ai;(r) = -27r(5(r) (2.1) 

Explicitly, v{r) = — ln(|r|/ro) where tq is a length scale set for simplicity to unity. We will 
work in the grand canonical ensemble, characterized by the inverse temperature [3 and 
the equivalent fugacities of the positively and negatively charged particles, z+ = Z- = z. 
Due to the charge ± symmetry, the induced particle densities are n_|_ = n_ = n/2 (n is 
the total number density of particles). 

Using the fact that — A/(27r) is the inverse operator of v{r) [see equation ( |2.1|) ], the 
grand partition function H of the 2D TCP can be turned via the Hubbard-Stratonovitch 
transformation (see, e.g., review [133]) into 



^ [Vl^p{-A{z))_ 
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where 



A{z) = J dV 



^ -2zcos{b(f)) 



.16.-" J ^'-'^^ 

6^ = (3/4 (2.3b) 

is the Euchdean action of the 2D classical sine-Gordon theory. Here, 0(r) is a real scalar 
field, / denotes the functional integration over this field and the fugacity z is renor- 
malized by a self-energy term. The sine-Gordon representation of the density of particles 
of one sign is 



j 

= z^ie'"'^) (2.4) 
where (• ■ ■) denotes the averaging over the sine- Gordon action ( p.3| ). The equality of 



densities n+ = for the considered z+ = Z- = z is a special case of a general symmetry 
relation (e^'^''^) = (e"'"*^), a arbitrary, which results from the invariance of the sine-Gordon 
action ( |2.3| ) with respect to the transformation —0. For two-body densities, one gets 

= ;2<,2,.(e''^'"^('-)e''^'''^('')) (2.5) 

Clearly, n++(r, r^) = n (r, r') and r2+_(r, r') = n |_(r, r^). For our purpose, it is useful 

to consider at the two-particle level also the pair distribution functions 

9aa'{r, r') = (2.6) 

the pair correlation functions 

K^ir,r')=g^^ir,r')-1 (2.7) 

and the Ursell functions 

U^a'ir, r') = n^n„,h^„>{r, r') (2.8) 
(denoted as p^^^T in refs. [|15|, [l^]). 



The parameter z in ( |2.3|) , i.e., the fugacity renormahzed by a diverging self-energy 
term, gets a precise meaning when one fixes the normahzation of the cos-field. In partic- 
ular, it was shown for the 2D TCP [^, that the behavior of the two-body densities for 
oppositely charged particles is dominated at short distance by the Boltzmann factor of 
the Coulomb potential, 

n+_(r, r') ~ 2;+2;_|r — r l^'^ as|r— r|^0 (2.9) 



With respect to (2.5) and the definition (|2.3bD , one has in the sine-Gordon picture 



^gi6<^(r)g-i6<^(r')^ ~ |r- r^]"^''' as |r- ^ (2.10) 

This short-distance normalization of the exponential fields makes the mapping complete. 
In quantum field theory, ( p. 10 ) is known as the conformal normalization, and it corre- 



sponds to taking (|2.3|) as a Gaussian conformal field theory perturbed by the relevant 
operator cos(60). The conformal normalization makes the bridge between the 2D TCP 
and the sine-Gordon model, and allows one to transfer results from one theory to the 
other |[. 

Let us now summarize the known facts about the sine-Gordon theory (|2.3|) . The theory 
has a discrete symmetry (l) + 2'Kj/h (j any integer). This symmetry is spontaneously 
broken, and the sine-Gordon theory is massive, in the region < 6^ < i (i.e., < /3 < 4 
for the 2D TCP). In this region, one has to consider one of infinitely many ground states 
{|0j)} characterized by = 27rj/6, say that with j = 0. The 2D classical sine-Gordon 
theory is integrable [^. Its particle spectrum consists of one soliton-antisoliton pair 
{A+,A_) (the corresponding particle index e = ± denotes the topological U{1) charge) 
with equal soliton and antisoliton masses m+ = m_ = M and of soliton-antisoliton bound 
states ("breathers") {Bfj = 1,2,... < 1/^} (the corresponding particle index e = j), 
whose number at a given 6^ depends on the inverse of the parameter 

The range < < 1 (0 < /3 < 2) is an attractive one (breathers exist), the range 
l<^<oo(2</3<4)isa repulsive one (breathers do not exist); ^ = 1 (/? = 2) 

6 



corresponds to the free-fermion (collapse) point. The masses of breathers {Bj} are given 
by the formula 

m^- = 2Msin l^yj^ (2.12) 

The lightest breather Bi is usually called the elementary one. The energy E and the 
momentum p of a particle e of mass are parametrized by the rapidity 6 G (—00, 00) as 
follows 

E = m^ cosh 9, p = m^smh9 (2.13) 

Like in other 2D integrable field theories, the A^-particle scattering amplitudes are purely 
elastic (the number of particles is conserved, the incoming and outgoing momenta are 
the same) and factorized into A^(A^ — l)/2 two-particle ^-matrices which are determined 
exactly by exploring their general properties: unitarity and crossing symmetry, validity 
of the Yang-Baxter equation and the assumption of "maximal analyticity" . 

The dimensionless specific grand potential ut of the sine-Gordon theory (^.2|), (p^.3|), 
defined by 

-co = lim — InS (2.14) 

y^oo V 

was found in ref. Q by using the Thermodynamic Bethe ansatz: 

mi 

= 7-7^ 2.15 

8sin(7rO 



Under the conformal normalization of the exponential fields ( |2.10|) , the relationship be- 



tween the parameter z and the soliton mass M was established in ref. 



r(6- 



,2^ 



7rr(l - 62) 



• v/?r((i + o/2) 
2r(e/2) 



H 2-262 



:2.i6) 



where F stands for the Gamma function. As a result, 



^ ' 2 dz 8z{l -h^) \2 j ^ ' 



Relations ( |2.14|) - ( p.l7|) , together with the equality 

n = 2z{e'^^) (2.18) 

constitute the complete set of equations for determining the thermodynamics of the 2D 
TCP Pj. They will be extensively used in the following section. 



3 Form-factor theory for the 2D TCP 

In 2D integrable systems, correlation functions of local operators Oa {a is a free parameter) 
can be written as an infinite convergent series over multi-particle intermediate states. For 
the two-point correlation function, one has 



N 



X 



^iF,,(^^,...,^i)exp -|r-r^|E"^.,cosh% (3.1) 



The first = term of the series is nothing but the decoupling {Oa){Oa')- The form 
factors 

i^a(^l,...,Mei...e^ = (0 1 a( 0) | (^1 ),•.•, (M ) (3.2a) 

^-■^^F,,(^^,...,^i) = (Z,^(M,...,Z,,(^^i)|a'(O)|0) (3.2b) 

are the matrix elements of the operator at the origin, between an A^-particle in-state 
(being a linear superposition of free one-particle states \Z^{9))) and the vacuum. They 
depend on particle rapidities only through their differences 6jk = 6j — 6k- The form 



factors can be determined exactly (for the underlying sine-Gordon theory, see refs. [p^] - 



30(1). They satisfy a set of constraint functional equations, originating from general 
properties of unitarity, analyticity, relativistic invariance and locality. The important fact 
is that these functional equations do not refer to a specific operator O, the S'-matrix is 
the only dynamical information needed. The nature of the operator is reflected only via 
normalization constants for the form factors. In what follows, under Oa we will understand 
an exponential operator 

a(r) = exp(ia0(r)) (3.3) 

The form-factor representation of two-point correlation functions ( p.l|) is particularly 
useful for large distance \r — r^\. In the limit |r — oo, the dominant contribu- 

tion to the truncated function (Ca('")C^a'('^))T = {Oa{r)Oa'{r')) — {Oa){Oa') comes from 
a multi-particle intermediate state with the minimum value of the total particle mass 
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Y.j=i f^nj) at the point of vanishing rapidities 9j 0. The corresponding exponential de- 
cay exp(— |r— r^ \ YljLi ^ej) is modified by a slower (inverse power law) decaying function 
which particular form depends on the form factors. For the sine-Gordon theory ( |2.3| ) (or, 
equivalently, the 2D TCP), in the region 0<6^<l/2(0</9<2), the lightest particle in 
the spectrum is the elementary breather Bi. As is clear from ( p.l2|) , approaching —* 1/2 
resp. ^ 1 (/3 ^ 2), its mass mi — > 2M. Note that, according to eq. ( p. 161) , M is finite 
at b"^ = 1/2, resp. ^ = 1 for a fixed value of z. The elementary breather Bi disappears at 
the free-fermion point b"^ = 1/2 (collapse point (3 = 2). The only existing particles remain 
the topologically neutral soliton-antisoliton pair {A^.,A_), which has just the same total 
mass 2M. Thus, the total mass remains 2M (which is temperature-dependent) in the 
whole region 1/2 < 6^ < i (2 < /3 < 4). 

3.1 < 6^ < 1/2 (0 < < 2) 

For the elementary breather _Bi, the form factors FaiOi, . . . , 6n)i...i and ^■■■^Fa'iON, . . . ,9i) 



= Fa'{0]y, . . . , ^^i)i...i are presented for the exponential operator (|3.3| ) in ref. In the 

notation (|3.2|), they read 

{0\e"''^\B,{9)) = -i {e"''^) (TrA)^/^ ''^^''(t'/)^^ (3.4a) 



(O|e^'^'^|i?i(02),Si(^i)) = -(e''^'^) (vrA) 

etc., where 



sin(7r^a/6)'' ^ 
sin(7r^) 



A = ^ sin«) cos ^ exp I - ^ 



dt t 



Ri9i - 02) (3.4b) 



(3.5) 



TT sin t J 

and the function R{6) is represented in the range — 27r + 7r^ < Im(^) < — tt^ by the integral 



R{e) 



TVexpjs^ 



exp < 4 



°° dt sinh t sinh(t^) sinh(t(l + ^)) 
~t sinh2(2t) 
dt sinh t sinh(t^) sinh(t(l +0)1 
~t sinh^(2t) J 



sinh^ t 1 



i^^ 



7r 



(3.6a) 
(3.6b) 
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An analytic continuation of this representation to the case of interest lm{6) = is given 
by the relation 

R{e)R{e ± Itt) 



sinh 6 



sinh 6 =1= i sin('7r^) 
Note a different notation in comparison with ref. |]5D| . 

The pair correlation function h(j„>{r, r') ( p.7|) is expressible as follows 



(3.7) 



(e 



icrb(f>{r) icr' b4>{r') 



(3i 



Using the form-factor representation ( |3.1| ), the dominant contribution to h„„r{r, r') in the 
limit \r— ^ oo is given by the one-breather Bi state with mass mi and the form factor 
( p.4a| ). Then, since 



d^ 



-rmicoshe _ Ko{mir) ~ ( 

\2mirJ 



1/2 



exp(— mirj 



(3.9) 



where Kq is the modified Bessel function of second kind [BH] , one finds, at asymptotically 
large distance, that 

haa'ij') ^ o-a'h{r) as r — > oo (3.10) 

with 

h{r) = —A exp(— mir) (3-11) 

\2mir J 

Using the thermodynamic formulae derived in section 2, the lightest-breather mass mi is 
expressed as follows 



mi 



sin(7r/3/(4-/3)) ^^/^ 
l-—(3'-—(3' + 0{P^) 



where 



K = {2Til3nf'^ 

denotes the inverse Debye length. The parameter A ( |3.5| ) takes the form 

7r/3 



A 



4 

— sm 

TT 



1 - 





1 cos 












384 



.2(4-/5) 



exp 



^ dt t 



TT sint 



96 2304 y ^ ^ 



(3.12a) 
(3.12b) 

(3.13) 

(3.14a) 
(3.14b) 
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The specific dependence of haa'{r) on tlie product of charges aa', formula (|3.10|) , means 
that the two-particle correlations are determined at large distances by the charge-charge 
correlation function. Indeed, 



J aa'K„,{r) = h{r) 



(3.15) 



On the other hand, the density-density correlation function Y.a,u'=± hacr'ij')/^ vanishes for 
the lowest one-breather Bi state, and becomes nonzero only for the next two-breather Bi 
state, with the form factor (|3.4b|) and much faster exponential decay exp(— 2mir). 



3.2 



1/2 < 62 < 1 (2 < < 4) 



The form factors for the soliton-antisoliton pair with topological f/(l)-charges e = ± 



are nonvanishing only for f/(l) neutral states Y^^=i^j = 0, with being inevitably an 



even number. Simultaneously, ''^•••^^Fa'(6'Ar, . . . ,9i) = Fa'{9N, ■ ■ ■ ,0i 



) ^lje„...e. 



where e* = — e. 



Form factors with different assignments of the charges {ej} are related by the relation 



Fa{Oi, . . . ,9i, 9i+i, . . . Om)^^. 



.eiei+l...ejv 



—Fa{9\, . . . , 6'i+i, 9i, . . . On)^^.. 



(3.16) 



and the hermiticity of the operator (|3.3|) implies 



Fa{6i, . . . , 6'Ar)ei...e^ = F^ai^N, • • • , 6'i)e^...6* (3-17) 



In the notation (|3.2| ), the two-particle form factors for the exponential operator ( |3.3| ) with 
the special value of a = 6 read |^ 



(3.18) 



where 



G{e) 



-2G{e) sinh(e)cotg exp 



^ sinh 



'6 + Itt' 



n -1 



exp 



t sinh(2t) cosht sinh(t4) \ tt^ 



(3.19a) 
(3.19b) 



Note a different notation. The form factors of the operator exp(— 160) are deducible from 
( CT ) and (|3l9D by using relation ( ^JTl) . 
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In the region 2 < /? < 4, for fixed the total particle density n is infinite, so that only 
the fugacity z can be chosen as a legitimate parameter. Instead of h, which vanishes in 
this region, the Ursell function (1^) is considered. Using ( ^.ltj| )-( p.iy| ) in the form-factor 
representation ( p.l|) given by the contribution of the soliton-antisoliton pair, we get, as 
previously. 



as r ^ oo 



Uaa'{r) ~ aa'U{r) 
where U{r) is the leading asymptotic r — oo term of the integral 
M4cos(7r/0 d0id^2 



(3.20) 



1664 



-G{e)G{-e) sinh'e 



X 



smh smh 



-Af r (cosh ei +cosh 6*2 ) 



(3.21) 



with e = 61-62. 

At /3 = 2, since ^ = 1 it holds G{6) = 1, and therefore 

.3 



U{r) 



m 



inr 



exp(— 2mr) 



P = 2 



where 



m = 2tiz 



(3.22a) 



(3.22b) 



is the soliton mass M at /3 = 2, see formula ( p.l6|) with 6^ = 1/2 and ^ = 1. 
For /? > 2, it can be readily shown from ( |3.21D that 



U{r) oc — rexp(-2Mr) 



2 < /? < 4 



(3.23a) 



where the soliton mass M is given by formula (|2.16|) in the following form 



M 



2r(/3/[2(4 -/?)]) 



'tizT{1- (3/A) 

r(/3/4) 



2-/3/2 



(3.23b) 



v/^r(2/(4-/3)) 

The prefactor in ( p.23a| ) could, in principle, be obtained from ( |3.21|) , but we omit this 
complicated calculation. 



12 



4 Analytic checks of the results 



4.1 Small coupling expansion 



In this section, the small-/? expansion ( pj.l2D and ( |3.14D of the asymptotic form ( |3.11D of the 



charge correlation function will be checked by a direct calculation using the renormalized 
Mayer expansion ||, . 

The two-dimensional Fourier transforms which will be needed are defined here, for 
instance for the charge correlation function h{r), as 



oo 



h(k) = / h(r)jQ(kr)rdr (4.1a) 
Jo 

/■oo „ 

h{r) = / h{k)Jo{kr)kdk (4.1b) 

JO 

where Jq is a Bessel function. h{k) is related to the Fourier transform c{k) of the charge 
direct correlation function by the Ornstein-Zernike relation 

h(k) = ^^^L (4.2) 



where n is the total number density of the particles (in ref. [r^|, h and c were called h' 
and c', respectively). 

At the lowest order in 13 (Debye-Hiickel approximation), c(r) = /31nr, c{k) = —j3/k'^, 
and 

m = (4.3) 

where n is the inverse Debye length defined in eq. ( p.l3| ). Thus, 

h{r) = -PKoinr) (4.4) 



where Kq is a modified Bessel function. The asymptotic form of h{r) is p6 

h(r) = -(3 exp(-/tr) (4.5) 

m ag reement with (|3TT|) - ( glil) . 
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Further terms in the /5-expansion of c(r) are given by the renormahzed Mayer expan- 
sion in density. Graphs of the excess Helmholtz free energy contributing to c(r) are the 
ones which have their two root-vertices with on odd bond-coordination and their field 
vertices with an even bond-coordination. Up to order P'^, for a fixed value of k,, 

c{r)=(3\nr + (3^C3{r)+(3%{r) (4.6) 

where /3^C3 and /S'^c^ are defined by the diagrams [where a wavy line represents the renor- 
mahzed bond —PKo^Kr)]: 



P^csir) = (^^^ = -(LK^(^Kr) (4.7a) 

and 



/3%(r)= ^/J^ =-(^Ko{nr)nK-^ Jd\Kr')K^{Kr')K^{n\r-r^\) (4.7b) 

These diagrams are of order and respectively [this is apparent for C4(r) when the 

factor nK~'^ in front of the integral is replaced by l/(27r/5)]. c(r) has no term of order (3'^. 
The Fourier transforms of Cs{r) and C4(r) will now be shown to have convenient integral 
representations. From now on, k^^ will be taken as the unit of length. 
In the Fourier transform 

1 r°° 

^3(^) = -a / K'oir)Mkr)rdr (4.8) 

O JO 

the integral can be viewed as l/(27r) times the scalar product of K^i^r) and Ko{r)Jo{kr), 
which can be written as the scalar product of their Fourier transforms. These transforms 
are [Sg 

In 




oo 



G{1) = I K^{r)Jo{lr)rdr 

and 



" ^./. , fl^^ 



(4.9) 



Hkil)= Ko(r)Jo(kr)Jo(lr)rdr = + + 1^ -2k^l^ + 2P + 21^)-^/^ (4.10) 
Jo 
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Therefore, setting / = 2sinh<^, one finds 
csik) = -i r G{l)Hk{l)m 

JO 

^3io 



dip if 



1/2 



(4.11) 



16 sinh^ + 16 sinh^ - 8(A;2 + 1) sinh^ + (A;2 + 1)2 

In the expression ([4.7b| ) of C4(r), there is an integral /(r), which is the convolution of Kq 
with itself. Its Fourier transform is /(/) = 27r[G'(/)]^. Again, in 



1 

Ci{k) = —— I Ko{r)I{r)Jo{kr)rdr 
Svr Jo 



(4.12) 



the integral can be viewed as 1/ (27r) times a scalar product, now of /(r) and KQ{r)Jo{kr) 
As above, in terms of the Fourier transforms /(/) and Hk{l)-, 



C4(fc) = -— / I{l)Hk{l)ldl = -- / [G'(/)]^iffe(/)M/ 

OTT JO 4 JO 



(4.13) 



Using (|4.9| ) and ( 4.10 ) in ( [4.13 ), and again setting I = 2sinhy9, one finds 

d(f (f"^ 



Ci{k) 



1 

'a Jo 



sinh ip cosh if 16 sinh^ (/? + 16 sinh^ y? — 8(^2 + 1) sinh^ ip + (A;2 + 1)2 
Finally, up to order 



1/2 

.14) 



c(fc) = --^[1 - /32A;%(A;) - /5=^A:2c4(fc)] 



(4.15) 



where 63 and C4 are given by ( [4.11|) and (|4.14|) . 
Using (|4.15|) in (|]2|) (with k = 1) gives 

f5[l - (5^k^c^{k) - p^eci{k)] 



h{k) 



(4.16) 



A;2 + l-/?2A;2c3(A;)-/?3pc4(A;) 

We are interested in the asymptotic behavior of h{r), which is governed by the poles of 
h{k) closest to the real axis. When j3 -—>■ 0, these poles are at A; = ±i, i.e., at k'^ = —1. For 
finding the location and the contributions of these poles in a /5 expansion when the 
and (3^ terms in the denominator of ( [4. 16] ) are taken into account, it is enough to expand 
C3 and C4 around /c^ = — 1 up to first order in 1 + fc^. 
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However, some care is needed for deriving this expansion from the expressions ( |4.11|) 
and ( [4.14|) . Indeed, these expressions were obtained by using the Fourier transform (|4.10|) 
of KQ{r)jQ{kr). However, for complex k, this integral diverges at large r when |Imfc| > 1. 
Thus, the integral representations ( [4 .111) and ( [4.14| ) are valid for 1 + k"^ real positive but 
are expected to exhibit some singularity at 1 + fc^ = [although the original function 
C3(/c) defined by ( [4. 81) and its counterpart for C4(/c) ( |4.12| ) are regular around 1 + fc^ = 0; 
for instance, ( [4. 81 ) defines a function of k which is analytical in the whole strip |Imfc| < 3, 
thus in particular for any real value of 1 + A;^ > —8] . The singularity in ( |4.11| ) or ( |4.14| ) 
can be seen directly if a naive expansion with respect to x = 1 + fc^ is attempted: the 
coefficient of is an integral which diverges at small (p. This is a warning that one 
must be careful when computing the coefficient of the previous term (of order x) in the 
expansion, as follows. 

The zeros of the denominator in ( [4.11|) are sinh^ = — (1/2)[1 — {x/2)] ± (1 — x)^/^], 
i.e., for small x, 

dip p 



(4.17) 



[sinh^ p + (xVl6) + . . .]i/2[cosh2 p - (x/2) + . . .]i/2 
While the second factor in the denominator can be expanded with respect to x, the first 
factor is the dangerous one in which the x^ term will give a contribution of first order in 
X and should not be suppressed. To first order in x. 



dp p 



12 Jo cosh(y9[sinh + (x^/lG)]-*^/^ 



1 + 



X 



4 cosh p 



dp sinh(y9 



12 70 cosh(/?[sinh> + (xVl6)]i/2 



+ 



(sinhy9 — p)dp 



12 Jo cosh (/?[sinh^ p + (xVl6)] V2 

X 



pdp 



(4.18) 



48 Jo cosh^ (/?[sinh^ p + (xVl6)]i/2 
In an expansion of C3 to first order in x, the third integral in the rhs of ( [4.18|) is needed 
only for x^ = 0. By an integration per partes, it is found to be 



dp p 



TV 



cosh p sinh p 8 



(4.19) 
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A naive expansion of the second integral in the rhs of (f4.18|) gives a finite coefficient for 
the term. Thus, the integral has no term of first order in x and can be evaluated at 



X 



0. Using tabulated integrals [BH], one finds 



dip (sinh (f — (f) 71 71^ 



(4.20) 



cosh {p sinh ip 

Finally, the first integral in the rhs of ( |4.18|) , which is the delicate one, can be exactly 
evaluated, by taking coshy) as the integration variable. The result does have a term of 
order x, which would have been missed if had been neglected in the integrand: 

dip sinhy? 



vr X 



cosh(y9[sinh + 



+ . . . 



(4.21) 



It should be noted that ( |4.21|) has been derived for x real positive. For x real negative, x 
must be replaced by |x| in the rhs of ([4.21|) , which confirms that the integral representation 
( [4.11| ) has a singularity (a kink) at x = 0. However, since the original expression ([4.8|) 



is regular around x = and the integral representation ( |4.11| ) is valid for x real positive, 
( [4.21| ) can be used for obtaining the correct analytical behavior of c^lk) near 1 + A;^ = 0. 
Using (|4l9| ), ( CT ), and (^?2l\) in ( ^ISD gives the 1 + fc^ expansion 



C3{k) 



96 



TT 



32 384 



(4.22) 



A similar approach can be used for expanding the expression (|4.14| ) of C4{k). To first 
order in x = 1 + A;^, 



dip sinh(y9 



+ 



16 Jo cosh^ (^[sinh^ Lp + (x2/16)]i/2 
1 r°° dip (sinh^ Lp - ip'^) 



16 Jq sinh^ ip cosh^ ip 
X 1"°° dip> 'pP' 



64 7o sinh^ ip cosh^ Lp 
By simple manipulations, the third integral in the rhs of ( |4.23| ) is found to be 



(4.23) 



dip) ipp- 



TT^ 1 



sinh y3cosh 36 3 



(4.24) 
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From tables the second integral is 

dip (sinh^ (f — (p^) 



sinh^ Lf cosh^ y9 



1-^ 
12 



(4.25) 



Finally, the first integral can be exactly evaluated and afterwards expanded in x with the 
result 

d(^ sinh ip , X , (^^26) 



— ^ 



/o cosh (y9[sinh + (a;^/16)] 
Using (|CT1 ), ( p5| ), and ( ]i:26| ) in ( ^:23D gives 

,2 



C4(/i;) 



TT 



1 



TT" 



(1 + A;2) + 



(4.27) 



192 V96 2304^ 

By using ( [4.22|) and ( |4.27|) in the denominator of ( |4.16| ), it is seen that, at order j3^, 
this denominator has poles aX k = ±imi with 



mi = 1 — 



192 



384 



(4.28) 



The asymptotic behavior of h{r) is governed by these poles and the corresponding residues, 
i.e., by writing h{k) in the form — A(A;^)(A;^ + m^)^^ and replacing k"^ by —m\ in X{k'^). 
One obtains 

MA:) ~ -^^^ (4.29) 

where, at order 



A = /3 



237rn 



^ 1^32 ^ 384 V96^ 2304j'^ 



(4.30) 



The corresponding h{r) is — A/Co(^i^) which has the asymptotic behavior (in units of 



h{r) 



TT 



2mi r 



1/2 



exp(— mirj 



(4.31) 



with uii and A given by ( [4.28|) and ([4.30 ). This is in agreement with the small-/? expansions 
of these quantities ( p.l2b| ) and (|3.14b| ) obtained from the exact result. 
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4.2 Free-fermion (3 = 2 point 



At the collapse point (3 = 2, the Ursell functions were found in refs. |T^, |17|. For fixed 
fugacity they read 

— j Klimr) (4.32a) 
— j Kl{mr) (4.32b) 
where m = 27iz, Kq and Ki are modified Bessel functions and a = ±. Consequently, 

Tfi 

Uaa'ir) ~ — acr' exp(— 2mr) as r — oo (4.33) 

Svrr 



in full agreement with ( p.20| ) and (|3.22|) . 



5 Conclusion 

In this paper, we have derived the exact large-distance behavior of particle correlation 
functions for the 2D two-component plasma by using the form-factor theory of the equiva- 
lent sine-Gordon model. The asymptotic decay is always exponential, with a continuously 
varying correlation length. The correlation length is determined by the particle spectrum 
of the sine-Gordon model: in the stability regime of the inverse temperatures < /3 < 2, 
it is equal to the inverse mass 1/mi of the lightest breather Bi, and in the collapse region 
2 < /3 < 4, where the breathers disappear, it is equal to the inverse total mass 1/(2M) of 
the only present particles, the soliton-antisoliton pair {A^,A_). On the other hand, the 
prefactor inverse-power-law function changes disco ntinuously at the collapse point /5 = 2: 
it behaves like l/i/r in the stability region < < 2 [see formula ( |3.11| )], like 1/r at 



/3 = 2 ( p.22| ) and like for 2 < /3 < 4 ( p.23|) . The change in the large-distance behavior 



of correlation functions from both sides of point /? = 2 is a sign of a singularity, which 
probably prevents the construction of an analytic (/3 — 2) expansion of Ursell functions and 
other finite statistical quantities, although all multi-particle Ursell functions are known 
at the collapse point . 
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The region 2 < /? < 4, in which the thermodynamic collapse occurs, deserves attention. 
The Ursell functions calculated in this region are equivalent to the corresponding Ursell 
functions of the lattice version of the 2D TCP, in the continuum limit when the lattice 



constant 0, as was done at /5 = 2 in refs. |jT5|, [T^. The phenomenon of the Kosterlitz- 



Thouless phase transition HTJ requires the fundamental presence of a maybe small, but 
nonzero, hard core attached to the charged particles. This is why we do not observe a 



typical K-T divergence of the correlation length oc exp(c/ as /9 4 . On the 
contrary, since the soliton mass M ( p.23b| ) goes to infinity as /3 — 4^ for fixed z, and even 



for z ~ 4 — /5 as was considered in the usual analysis of the renormalization flow close 
to the K-T transition, the correlation length goes to zero when approching (3 = 4. In a 
certain sense, likewise as /5 = 2 is the collapse point for the thermodynamics, (3 = 4 is the 
"collapse" point for the Ursell correlation functions of our model of pointlike particles. 
In the interval < (3 < 2, the introduction of a small hard core is a slight perturbation 
which does not change the thermodynamics and the correlation functions substantially. 

A next natural step is to study inhomogeneous situations when the 2D TCP is con- 
fined by an impermeable wall of dielectric constant ew- In such situations, the asymptotic 
decay of pair correlation functions is direction-dependent. In the case of ideal conductor 
{ew — * oo) and ideal dielectric {e^ = 0) rectilinear walls, the plasma can be formulated 
as a boundary sine-Gordon model with Dirichlet [0 and Neumann |^ boundary condi- 
tions, respectively. One of first attempts to determine the form factors of such integrable 
boundary theories was made in ref. . 
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